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I.  INTRODUCTION 


Target  configurations  for  studying  penetration  by  shaped-charge 
jets  are  becoming  more  complex;  consequently,  numerical  techniques, 
capable  of  performing  computations  with  more  than  two  materials  in  a 
computational  cell  (such  as  the  case  with  the  DORF1  and  TRIDORF2  codes) , 
are  required.  The  advanced  operator  splitting  versions  of  DORF  and 
TRIDORF,  referred  to  as  the  2DS0IL  and  3DS0IL  codes,  respectively,  have 
the  capability  of  treating  up  to  four  materials  within  a  computational 
cell . 


An  accurate  treatment  of  free-surface  movement  is  necessary  for 
thin  targets  and  multiple  material  targets.  Special  free-surface 
treatments  have  been  incorporated  into  the  2DS0IL  aud  TRIDORF  codes. 

For  computational  problems  involving  large  stresses,  the  available 
viscosity  due  to  material  advection  is  adequate  to  accurately  describe 
the  stress,  density,  and  velocity  fields.  However,  at  lower  stresses, 
an  artificial  viscosity  prescription  is  required  to  dampen  the 
oscillatory  flow  results.  An  artificial  viscosity  routine  which  is 
incorporated  into  the  solution  during  compression  and  rarefaction  phases 
is  now  available  in  the  2DS0IL  and  3DS0IL  codes. 

The  early  versions  of  TRIDORF  and  3DS0IL  contained  very  primitive 
and  limited  methods  for  generating  starting  data  for  the  codes.  Generator 
codes  (similar  to  the  CLAM2  code)  have  been  developed  and  checked  out 
for  the  3DS0IL  and  TRIDORF  codes. 

Test  problems  (Section  VI,  Summary  and  Conclusions)  have  been 
completed  for  the  purpose  of  comparing  two-dimensional  results  from  a 
two-dimensional  code  to  the  two-dimensional  results  from  a  three-dimensional 
code. 


Sections  II  thru  V  contain  operating  instructions  for  the  Ballistic 
Research  Laboratory1 s  versions  of  the  2DCLAM,  2DS0IL,  3DCLAM,  and  3DS0IL 
codes . 


1  Johnson,  W.  E.,  "Development  and  Applications  of  Computer  Programs  to 
Hypervelocity  Impact ,"  Systems,  Science  and  Software  Report  3SR-749, 

July  1971. 

2 Johnson,  W.  E.,  "TRIDORF  -  A  Two-Material  Version  of  the  TRIOIL  Code 
with  Strength,"  Computer  Code  Consultants  Report  CCC-976,  September  1976. 
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II.  2 DC LAM  INPUT 


Units  are  in  jerks,  grams,  kev,  centimeters,  and  shakes  where 
1  jerk  =  1016  ergs  and  1  shake  =  10~8  s. 


An  asterisk  signifies  that  the  datum  is  entered  in  a  real  number 


format;  otherwise, 

it  is  entered  in 

an  integer  format. 

CARD  NUMBER 

COLUMN 

DESCRIPTION 

1 

2-72 

Header  card  for  comments 

2 

*  1-10 

Problem  number 

*  11  -  20 

Maximum  number  of  cells  in  the 
r-direction  (IMAX) 

*  21  -  30 

Maximum  number  of  cells  in  the 
z-direction  (JMAX) 

*  31  -  40 

0.  for  hydrodynamic  and  radiation 
computations  or 

1.  for  hydrodynamic  and  strength 
computations 

*  41  -  50 

0.  for  no  tracer  particles  or 

2.  for  tracer  particles 

*  51  -  60 

0.  for  r-z  geometry  or 

1.  for  x-y  plane-strain  geometry 

*61-70  Calls  the  REZONE  routine  and  is 

equal  to  the  tape  number  from  ADJUST; 
otherwise,  set  to  0. 

71  -  72  Binary  tape  number 


3  (Cell  Dimensions) 
Several  Number  type  cards 
are  required. 


1  1  if  this  is  the  last  DX  or  DY  card  or 
0  if  there  are  more  DX  or  DY  cards 

2  0  for  DX  data  or 
1  for  DY  data 


3  -  4  Number  of  cells  that  will  have  the 

DX  or  DY  value  in  columns  11  -  20 


5  -  6  Number  of  cells  that  will  have  the 

DX  or  DY  value  in  columns  21-30 
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CARD  NUMBER 


COLUMN 


DESCRIPTION 


7-8 

9-10 

*  11  -  20 

*  21  -  30 

*  31  -  40 

*  41  -  50 

If  DX  or  DY  above  is  negative, 
DY (I)  -  DY(I  -  1)* *ABS (value). 

Last  card  of  set  *  21  -  30 

Insert  data  cards  for  CUBIT  if 

Geometry  and  p,  I,  U,  and  V  Input 

CARD  NUMBER  COLUMN 

1  1 

5-7 

*  11  -  20 

*  21  -  30 

*  31  -  40 

*  41  -  50 

*  51  -  60 


Number  of  cells  that  will  have  the 
DX  or  DY  value  in  columns  31  -  40 

Number  of  cells  that  will  have  the 
DX  or  DY  value  in  columns  41  -  50 

Value  of  DX  or  DY 

Value  of  DX  or  DY 

Value  of  DX  or  DY 

Value  of  DX  or  DY 

DX(I)  =  DX(I  -  1)*ABS (value)  or 

80. 

S9  >  0. 

DESCRIPTION 

Load  a  1 

N2,  the  number  of  particles  per  cell 
(1  <_N  <  20) 

Z -coordinate  for  the  origin  of  the 
radius  vector  for  the  FIT  routines 
(YC) 

R-coordinate  for  the  origin  of  the 
radius  vector  for  the  FIT  routines 
(XC) 

FIT  number  (1  thru  6)  to  use  to 
compute  I,  U,  V,  and  p 

Material  number  (1  thru  4) 

Specific  heat 
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Following  the  first  card  of  each  package  are  five  possible  data 
cards : 

1 .  Generate  geometry 

2.  Non-generate  geometry 

3.  Density 

4.  Specific  internal  energy 

5.  Velocity  (U  and  V) 


Note:  Many  Number  1  and  Number  2  cards  may  be  used  (maximum  of  72), 
but  only  one  each  of  cards  3,  4,  and  5  may  be  used. 


GEOMETRY 


1.  Rectangle 


2.  Triangle 


3.  Ellipse  or  Circle 


COLUMN 

1 

7 

*  11  -  20 

*  21  -  30 

*  31  -  40 

*  41  -  50 

1 

7 

*11-20 

*21-30 

*31-40 

*  41  -  50 
*51-60 

*  61  -  70 

1 

2 
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DESCRIPTION 


Load  a  4 


Load  a  1  to  generate 
Load  a  0  to  void 

Left  X  of  rectangle 

Right  X  of  rectangle 

Lower  Y  of  rectangle 

Upper  Y  of  rectangle 

Load  a  6 

Load  a  1  to  generate 
Load  a  0  to  void 

Xl-coordinate  of  vertex  1 

Yl-coordinate  of  vertex  1 

X2-coordinate  of  vertex  2 

Y2-coordinate  of  vertex  2 

X3- coordinate  of  vertex  3 

Y3-coordinate  of  vertex  3 

Load  a  4 

Load  a  1 

Load  a  1  to  generate 
Load  a  0  to  void 
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GEOMETRY 


COLUMN 


DESCRIPTION 


*11-20  Semi-axis  in  R  of  an  ellipse; 
radius  of  a  circle 


*21-30  Semi-axis  in  Z  of  an  ellipse; 
0.  for  a  circle 


*31-40  R  coordinate  of  the  center  of  the 
ellipse  or  circle 


*41-50  Z  coordinate  of  the  center  of  the 
ellipse  or  circle 


Following  the  geometry  cards  are  data  cards  that  refer  to  all  cells 
within  the  package : 


1.  Density  card  (load  a  51  in  columns  1-2) 

2.  Specific  internal  energy  card  (load  a-52  in  columns  1-2) 

3.  Velocity  (U  and  V)  card  (load  a  53  in  columns  1-2) 


The  data  for  p  (TABR  array) ,  I  (TABI  array) ,  U  (TABUV  array) ,  and 
V  (TABUV  array)  are  loaded  in  real  number  format  in  columns  11  -  20, 

21  -  30,  31  -  40,  41  -  50,  51  -  60,  and  61  -  70. 

FIT  routines  are  provided  for  analytical  or  table  fits  for  each 
package.  FIT3,  FIT4,  FIT5,  and  FIT6  are  blank  and  are  available  to  the 
user.  FIT1  and  FIT2  are  described  below.  FIT  numbers  are  specified  on 
the  first  card  of  each  package. 

The  FIT  routine  will  assign  a  density,  specific  internal  energy, 
and  velocity  components  to  each  particle,  N,  having  coordinate  TTX  (the 
r-coordinate  of  particle  N  relative  to  the  origin  of  the  radius  vector) 
and4 coordinate  TTY  (the  z-coordinate  of  particle  N  relative  to  the  origin 
of  the  radius  vector) . 

The  standard  output  from  the  FIT  routines  is 


WSR,  the  particle  density  of  particle  N 

WSI,  the  specific  internal  energy  of  particle  N 

WSU,  the  r- velocity  component  of  particle  N 

WSV,  the  z-velocity  component  of  particle  N 


FIT1 


FORTRAN  Expression 


Analytic  Expression 


WSR  =  TABR(l) 
WSI  =  TABI ( 1) 
WSU  =  TABUV (1) 
WSV  =  TABUV(2) 


p  =  constant 
I  =  constant 
U  =  constant 
V  =  constant 
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FIT2 


FORTRAN  Expression 

WS  =  SQRT (TTX**2  +  TTY**2) 
WSR  =  TABR(l)  +  TABR (2) *WS 
WSI  =  TABI(l)  +  TABI (2) *WS 
WSS  =  TABUV(l)  +  TABUV (2) *WS 
WSU  =  WSS*TTX/WS 
WSV  =  WSS*TTY/WS 


Analytic  Expression 

R  =  (r2  +  z2)^ 
p  =  a  +  bR 
I  =  c  +  dR 
S  =  e  +  fR 
U  =  Sr/R 
V  =  Sz/R 

where  a,  b,  c,  d,  e,  and  f  are 
constants . 


The  last  card  of  the  geometry  set  will  have  a  2  in  column  1. 
Additional  data  cards  will  be  called  if  N10  >  1. 


Tracer  Particle  Option 

COLUMN  1  COLUMNS  4-6  VARIABLE  DESCRIPTION 


2 

56 

LI 

The  number  of  particle  packages  to 
be  generated 

1 

100 

14 

Load  a  1.0 

There 

will  be 

LI  of  the 

following  sets: 

2 

67 

N1 

Number  of  tracer  particles  for  thi: 
set  (for  a  circle,  the  code  gives 
2*N1  particles 

1.  Circle 

Option 

1 

29 

SI 

0.  for  a  circle  of  particles 

1 

30 

S2 

R-coordinate  of  the  center  of  the 
circle 

1 

31 

S3 

Z-coordinate  of  the  center  of  the 
circle 

1 

32 

S4 

Radius  of  the  circle 

2.  Line  Segment  Option 

1 

29 

SI 

1.  for  a  line  segment  of  particles 

1 

30 

S2 

Xl-coordinate  of  an  end  point  of 
the  line 
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COLUMN  1  COLUMNS  4-6  VARIABLE  DESCRIPTION 


1 

31 

S3 

Yl-coordinate  of 
the  line  segment 

an 

end  point 

of 

1 

32 

S4 

X2-coordinate  of 
the  line  segment 

an 

end  point 

of 

1 

33 

S5 

Y2-coordinate  of 
the  line  segment 

an 

end  point 

of 

Each  of  the  LI  sets  will  end  with  a  card: 


COLUMN  1 

COLUMNS  4-6 

COLUMNS  8-9 

1 

100 

1. 
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III.  2DS0IL  INPUT 


Units  are  in  jerks,  grams,  kev,  centimeters,  and  shakes  where 

1  jerk  =  1016  ergs  and  1  shake  =  10"8  s. 

All  data  loaded  via  the  CARD  routine  are  in  a  real  number  format. 
A  2  in  column  1  instructs  the  CARD  routine  (data  loader)  to  convert 
input  real  numbers  to  an  integer  format. 

There  are  three  sets  of  data  cards.  A  1  in  column  1  designates  a 
card  to  be  the  last  in  the  second  and  third  sets  of  data  cards. 

First  Set 


Card  1:  Column  3  contains  N,  the  number  of  header  cards  that  are 
used  for  comments. 

Next  N  cards:  Comments  in  columns  2  thru  72. 


Second  Set 


Column  1 

Columns  4-6 

Variable 

Description 

2 

73 

N7 

Number  of  binary  tape  to  be  read 

2 

85 

NMAX 

0.  for  no  tracer  particles  or 

1.  for  tracer  particles 

2 

96 

NST 

0.  for  hydrodynamic  computations  or 
1.  for  hydrodynamic  and  strength 
computations 

151 

PK(1) 

Problem  number 

1 

152 

PK(2) 

Cycle  number  for  restart 

Third  Set 

Column  1 

Columns  4-6 

Variable 

Description 

6 

PMIN 

Minimum  pressure 

9 

UMIN 

Minimum  velocity 

10 

EIMIN 

Minimum  specific  internal  energy 

12 


Column  1 


Columns  4-6 


Variable 


Description 


11 

RADC 

0,  for  no  radiation  computations  or 
* 0137  jerks/cc/kev  for  radiation 
diffusion  computations 

16 

THEMIN 

Minimum  temperature  (kev) 

17 

CLITE 

300  cm/shake,  radiation  velocity 

18 

ALPHA 

1.,  implicit  coefficient  for 
radiation  flow 

21 

DTOPTN 

-1*,  the  code  will  control  the  time 
step  using  the  Courant  and  particle 
velocity  conditions  or 

1.,  the  code  will  use  the  loaded 
constant  time  step 

23 

ESTAB 

Early  time  stability 

24 

EPSLN 

Minimum  flux  (p)  at  free  surface 
when  E  >  EMELT 

25 

FSTAB 

Final  time  stability 

26 

RHOMIN 

Minimum  density  for  cell  to  be 
checked  for  time  change  stability 

27 

BGTHE 

Background  temperature  (kev) 

28 

TCUT 

Use  units  of  04 

If  9  >  TCUT,  the  stability  (radiation) 
is  checked  if  the  zone  is  either 
cooling  or  heating,  or, 
if  0  < TCUT,  the  stability  (radiation) 
is  checked  only  if  the  zone  is  cooling 

29 

SI 

.5,  factor  in  MFLUX  routine 

37 

S9 

.01,  epsilonics  in  strength 

38 

S10 

Projectile  velocity  at  t  =  0. 

43 

FEF 

allowable  energy  change,  AE/E,  per 
cycle  for  radiation  flow 

44 

FACT 

Factor  for  increasing  ESTAB  per  cycle 
until  ESTAB  =  FSTAB 

13 


Column  1 


Columns  4-6 


Variable 


Description 


46 

DMIN 

Problem  will  stop  on  energy  check 
if  | AE/E/cycle |  <  DMIN 

47 

Z(47) 

Constant  used  in  flux  limiter  at  free 
surface  in  the  radiation  flow 

49 

2(49) 

Minimum  DX  or  DY 

Note:  LI  and  L2  are  right  boundary  values  of  a  cell.  L3  and  L5  are  top 
boundary  values  of  a  cell.  L4  is  the  bottom  boundary  value  of  a  cell. 


2  54 


2  55 


2  56 


2  57 


2  58 


LI  Active  grid  counter  that  is  advanced 
or  decreased  after  a  computational 
cycle 

L2  Active  grid  counter  that  is  advanced 
or  decreased  after  a  computational 
cycle 

L3  Active  grid  counter  that  is  advanced 
or  decreased  after  a  computational 
cycle 

L4  Active  grid  counter  that  is  advanced 
or  decreased  after  a  computational 
cycle 

L5  Active  grid  counter  that  is  advanced 
or  decreased  after  a  computational 
cycle 


Note:  1.  indicates  a  reflective  boundary  condition; 

-1.  indicates  a  transmittive  boundary  condition. 


2 

59 

LHL 

2 

60 

LHT 

Left  hydrodynamic  boundary  condition 

Top  hydrodynamic  boundary  condition 
14 


Column  1 

Columns 

4-6  Variable 

Description 

2 

61 

LHR 

Right  hydrodynamic  boundary  condition 

2 

62 

LHB 

Bottom  hydrodynamic  boundary  condition 

2 

63 

LRL 

Left  radiation  boundary  condition 

2 

64 

LRT 

Top  radiation  boundary  condition 

2 

65 

LRR 

Right  radiation  boundary  condition 

2 

66 

LRB 

Bottom  radiation  boundary  condition 

2 

67 

NI 

Maximum  number  of  iterations  in  a 
mixed  cell  to  achieve  pressure 
equilibrium 

2 

75 

N9 

0.  for  no  strength  computations  or 

1.  for  strength  computations 

2 

86 

IFRAD 

Maximum  number  of  iterations 
(Newton-Raphson)  for  calculating  the 
temperature  from  the  energy  equation 

2 

96 

NST 

0.  for  no  strength  computations  or 

1.  for  strength  computations 

110 

Z(110) 

1.1*THETA  (background  for  radiation- 
hydrodynamic  computations)  or 

1.1*UMIN  (for  hydrodynamic-strength 
computations) 

146 

TLP 

Time  frequency  for  long  print 

147 

DELP 

Time  interval  for  long  print 

148 

TDUMP 

Time  frequency  for  tape  dump 

149 

DELTD 

Time  interval  for  tape  dump 

150 

TSTOP 

Time  for  problem  to  terminate  normally 

2 

161 

IDEQST  (1) 

Equation  of  state  number  for 
material  1 

2 

162 

IDEQST(2) 

Equation  of  state  number  for 
material  2 

2 

163 

IDEQST (3) 

Equation  of  state  number  for 

material  3 
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Column  1 

Columns  4-6 

Variable 

Description 

2 

164 

IDEQST(4) 

Equation  of  state  number 
material  4 

2 

170 

IDEQST (10) 

Total  number  of  materials  used  in 
this  problem 

2 

171 

IDKAP(l) 

Opacity  routine  for  material  1 

2 

172 

IDKAP (2) 

Opacity  routine  for  material  2 

2 

173 

IDKAP (3) 

Opacity  routine  for  material  3 

2 

174 

IDKAP (4) 

Opacity  routine  for  material  4 

2 

180 

IDKAP (10) 

0.  for  no  artificial  viscosity 
computations  or 

1.  for  artificial  viscosity 
computations 

181 

OPCON(l) 

Linear  viscosity  coefficient 

182 

OPCON (2) 

Quadratic  viscosity  coefficient 

183 

OPCON (3) 

Maximum  change  in  velocity.  If  the 
change  is  greater  than  OPCON (3),  the 
artificial  viscosity  coefficient  is 
set  to  0. 

Equation  of 

State  Data 

M  =  25 (N  -  1)  where  N  is  the  material  number  (1  thru  4). 

EQSTC (M+l)  thru  EQSTC(M+10)  are  constants  for  Tillotson's  equation 
of  state.  EQSTC (M+ll)  thru  EQSTC (M+25)  are  constants  that  are  required 
for  any  form  of  an  equation  of  state. 


Column  1  Columns 

4-6  Variable 

Description 

431+M 

EQSTC (M+l) 

P8  (g/cc) 

432+M 

EQSTC (M+ 2) 

a 

433+M 

EQSTC (M+3) 

J0  (jerk/g) 

434+M 

EQSTC (M+ 4) 

b  (jerk) 

435+M 

EQSTC (M+ 5) 

A  (jerk/cc),  bulk  modulus 
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Column  1  Columns  4-6  Variable 
436+M  EQSTC (M+6) 

437+M  EQSTC (M+7) 

438+M  EQSTC (M+8) 

439+M  EQSTC (M+9) 

440+M  EQSTC (M+ 10) 

441+M  EQSTC (M+ 11) 

442+M  EQSTC (M+12) 

443+M  EQSTC (M+13) 

444 +M  EQSTC (M+ 14) 

445+M  EQSTC (M+ 15) 

446 *M  EQSTC (M+ 16) 

447+M  EQSTC (M+ 17) 

448+M  EQSTC (M+18) 

449+M  EQSTC (M+ 19) 

450+M  EQSTC (M+20) 

451+M  EQSTC (M+21) 

452+M  EQSTC (M+22) 

453+M  EQSTC (M+23) 

454+M  EQSTC (M+24) 

455+M  EQSTC (M+25) 


Description 
Ig  (jerk/ g) 

P  (jerk/g) 

a 

3 

B  (jerk/g) 

Not  used 
Not  used 
Not  used 
Not  used 
io-£*p0 

10-6*p0 

IMELT  (jerk/g) 

Y0  (jerk/g),  yield  strength  in  shear 

a  in  Y  =  (YQ  +  aP) 

y,  shear  modulus 

p/pQ  failure  criterion 

Maximum  pressure  in  Y  function  above 

Not  used 

Not  used 

Specific  heat  (jerk/g/kev) 
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IV.  3 DC LAM  INPUT 


Units  are  in  jerks,  grams,  kev,  centimeters,  and  shakes  where 
1  jerk  =  1016  ergs  and  1  shake  =  10"8  s. 

An  asterisk  signifies  that  the  datum  is  entered  in  a  real  number 
format;  otherwise,  it  is  entered  in  an  integer  format. 

CARD  NUMBER  COLUMN  DESCRIPTION 


1 

2 

-  72 

2 

*  1 

-  10 

*  11 

-  20 

*  21 

-  30 

*  31 

-  40 

*  51 

-  60 

71  -  72 

3  (Cell  Dimensions)  1 

At  least  3  Number  3  type 
cards  are  required. 

2 

3-4 

5-6 

7-8 


Header  card  for  comments 
Problem  number 

Maximum  number  of  cells  in  the 
x-direction  (IMAX) 

Maximum  number  of  cells  in  the 
y-direction  (JMAX) 

Maximum  number  of  cells  in  the 
z-direction  (KMAX) 

0.  for  hydrodynamic  computations  or 
1.  for  hydrodynamic  and  radiation 
computations  or 

-1.  for  hydrodynamic  and  strength 
computations 

Binary  tape  number 

1  if  this  is  the  last  DX,  DY,  or  DZ 
data  card  or 

0  if  there  are  additional  DX,  DY,  or 
DZ  data  cards 

0  for  DX  data  or 

1  for  DY  data  or 

2  for  DZ  data 

Number  of  cells  that  will  have  the 
DX,  DY,  or  DZ  value  in  columns  11-20 

Number  of  cells  that  will  have  the 
DX,  DY,  or  DZ  value  in  columns  21-30 

Number  of  cells  that  will  have  the 
DX,  DY,  or  DZ  value  in  columns  31-40 
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CARD  NUMBER 


COLUMN 


DESCRIPTION 


*  11 
*  21 


10 

Number  of  cells  that 
DX,  DY,  or  DZ  value 

will  have  the 
in  columns  41- 

20 

Value  of  DX, 

DY, 

or 

DZ 

30 

Value  of  DX, 

DY, 

or 

DZ 

40 

Value  of  DX, 

DY, 

or 

DZ 

50 

Value  of  DX, 

DY, 

or 

DZ 

31 


If  DX,  DY,  or  DZ  above  are  negative,  DX(I)  =  DX(I  -  1) *ABS (value) , 

DY Cl)  =  DY(I  -  1 )*ABS (value) ,  or  DZ(I)  =  DZ(I  -  1) *ABS (value) ,  respectively. 

Last  card  of  set  *  21  -  30  80. 

Geometry  and  p,  I,  U,  V,  and  W  Input 

Data  is  loaded  to  generate  the  following  packages: 

1.  Parallelepiped 

2.  Tetrahedron  which  is  triangular  in  the  x-y  plane 
and  constant  for  all  z  planes 

3.  Sphere 

4.  Right  circular  cylinder  with  its  axis  parallel 
to  the  y-axis 


CARD  NUMBER 


COLUMN 

1 

5-7 


DESCRIPTION 


Load  a  1 


11 


21 


20 


30 


*31-40 


N3,  the  number  of  particles  per  cell 
(1  <  N  £  10) 

X-coordinate  for  the  origin  of  the 
radius  vector  for  the  FIT  routines 
(XC) 

Y-coordinate  for  the  origin  of  the 
radius  vector  for  the  FIT  routines 
(YC) 

Z-coordinate  for  the  origin  of  the 
radius  vector  for  the  FIT  routines 
(ZC) 


19 


CARD  NUMBER 


COLUMN 


DESCRIPTION 


*  41  -  50  FIT  number  (1  thru  6)  used  to 

compute  p,  I,  U,  V,  and  W  for  the 
N  particles 

*  51  -  60  Material  number  (1  thru  4) 

*  61  -  70  Specific  heat  for  material  number 

1,  2,  3,  or  4 

Following  the  first  card  of  each  package  are  five  other  types  of 
data  cards : 

1.  Generate  geometry  with  options  below 

2.  Delete  geometry  with  options  below 

3.  Density  (only  one  per  package) 

4.  Specific  internal  energy  (only  one  per  package) 

5.  Velocities  U,  V,  or  W  (only  one  per  package) 

The  density,  specific  internal  energy  and  velocity  components  are 
generated  per  particle  in  the  FIT  routines. 


For  data  cards  1  and  2  above,  3DCLAM  has  the  options  for  generating 
or  deleting  the  following  geometries: 


GEOMETRY 

COLUMN 

DESCRIPTION 

1.  Parallelepiped 

1 

Load  a  4 

7 

Load  a  1  to  generate  or 

load  a  0  to  delete  or  omit  generating 

*  11 

-  20 

Left  x-coordinate  of  parallelepiped 
(XI) 

*  21 

-  30 

Right  x-coordinate  of  parallelepiped 
(X2) 

*  31 

-  40 

Bottom  y-coordinate  of  parallelepiped 
(Yl) 

*  41 

-  50 

Top  y-coordinate  of  parallelepiped 
(Y2) 

*  51 

-  60 

Back  z-coordinate  of  parallelepiped 
(Zl) 

*  61 

-  70 

Front  z-coordinate  of  parallelepiped 
(Z2) 
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COLUMN 
1 
7 

*  11  -  20 

*  21  -  30 

*  31  -  40 

*  41  -  50 

*  51  -  60 

*  61  -  70 

3.  Sphere  1 

2 

7 

*  11  -  20 

*  21  -  30 

*  31  -  40 

*  41  -  50 

*  51  -  60 

4.  Right  circular  cylinder  1 
with  axis  parallel  to  the 

y- axis 

2 

7 

* 

*  11  -  20 


DESCRIPTION 

Load  a  6 

Load  a  1  to  generate  or 

load  a  0  to  delete  or  omit  generating 

X-coordinate  of  vertex  1 

Y-coordinate  of  vertex  1 

X-coordinate  of  vertex  2 

Y-coordinate  of  vertex  2 

X-coordinate  of  vertex  3 

Y-coordinate  of  vertex  3 

Load  a  4 

Load  a  1 

Load  a  1  to  generate  or 

load  a  0  to  delete  or  omit  generating 

Radius  of  sphere 

0. 

X-coordinate  of  the  center  of  the 
sphere 

Y-coordinate  of  the  center  of  the 
sphere 

Z-coordinate  of  the  center  of  the 
sphere 

Load  a  6 

Load  a  1 

Load  a  1  to  generate  or 

load  a  0  to  delete  or  omit  generating 

Radius  of  cylinder 


GEOMETRY 
2.  Tetrahedron 
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GEOMETRY 


COLUMN 


DESCRIPTION 


*21-30  X-coordinate  of  the  center  of  the 
cylinder 

*  31  -  40  Z-coordinate  of  the  center  of  the 

cylinder 

*  41  -  50  Y-coordinate  of  the  bottom  of  the 

cylinder 

*51-60  Y-coordinate  of  the  top  of  the 
cylinder 


Following  the  geometry  cards  are  the  following  data  cards  that  refer 
to  all  the  particles  within  this  package: 

1.  Density  card  (load  a  51  in  columns  1-2) 

2.  Specific  internal  energy  card  (load  a  52  in  columns  1-2) 

3.  Velocity  card  (load  a  53  in  columns  1-2) 

If  the  density,  specific  internal  energy,  and  velocities  for  this 
package  are  the  same  as  the  previous  package,  a  51,  52,  or  53  card  is 
not  required. 


The  data  for  p  (TABR  array),  I  (TABI  array),  U  (TABUV  array), 

V  (TABUV  array) ,  and  W  (TABUV  array)  are  loaded  in  real  number  format  in 
columns  11  -  20,  21  -  30,  31  -  40,  41  -  50,  51  -  60,  and  61  -  70. 


FIT  routines  are  provided  for  analytical  or  table  fits  for  each 
package.  FIT3,  FIT4,  FIT5,  and  FIT6  are  blank  and  are  available  to  the 
user.  FIT1  and  FIT2  are  described  below.  FIT  numbers  are  specified  on 
the  first  card  of  each  package. 


The  FIT  routine  will  assign  a  density,  specific  internal  energy, 
and  velocity  components  to  each  particle,  N,  having  coordinate  TX  (the 
x-coordinate  of  particle  N) ,  coordinate  TY  (the  y-coordinate  of  particle 
N) ,  and  coordinate  TZ  (the  z-coordinate  of  particle  N) .  If  point 
(TX,TY,TZ)  is  other  than  point  (0,0,0),  the  coordinates  of  particle  N 
are  relocated  to  point  (TTX,TTY,TTZ)  by 

TTX  =  TX  -  XC 
TTY  =  TY  -  YC 
TTZ  =  TZ  -  ZC 


where  point  (XC,YC,ZC)  is  the  origin  of  the  position  vector. 

The  standard  output  from  the  FIT  routines  is 

WSR,  the  particle  density  of  particle  N 

WSI,  the  specific  internal  energy  of  particle  N 

WSU,  the  x- velocity  component  of  particle  N 

WSV,  the  y- velocity  component  of  particle  N 

WSW,  the  z-velocity  component  of  particle  N 
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r 


FITl 


FORTRAN  Expression 

Analytic  Expression 

WSR  =  TABR(l) 

WSI  =  TABI(l) 

WSU  =  TABUV(l) 

WSV  =  TABUV(2) 

WSW  =  TABUV(3) 

p  =  constant 

I  =  constant 

U  =  constant 

V  =  constant 

W  =  constant 

FIT2 

FORTRAN  Expression 

Analytic  Expression 

WS  =  SQRT(TTX**2  +  TTY**2  +  TTZ**2) 
WSR  =  TABR(l)  +  TABR(2) *WS 

WSI  =  TABI(l)  +  TAB I (2) *WS 

WSS  =  TABUV(l)  +  TABUV(2) *WS 

WSU  =  WSS*TTX/WS 

WSV  =  WSS*TTY/WS 

WSW  =  WSS*TTZ/WS 

R  =  (x2  +  y2  +  z2)"2 
p  =  a  +  bR 

I  =  c  +  dR 

S  =  e  +  fR 

U  =  Sx/R 

V  =  Sy/R 

W  =  Sz/R 

where  a,  b,  c,  d,  e 
constants . 

The  last  card  of  the  geometry  set  will  have  a  2  in  column 


are 

1. 
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V.  3DS0IL  INPUT 


Units  are  in  jerks,  grams,  kev,  centimeters,  and  shakes  where 

16  *  6 

1  jerk  =  10  ergs  and  1  shake  =10  s. 

All  data  loaded  via  the  CARD  routine  are  in  a  real  number  format. 
A  2  in  column  1  instructs  the  CARD  routine  (data  loader)  to  convert 
input  real  numbers  to  an  integer  format. 

There  are  three  sets  of  data  cards.  A  1  in  column  1  designates  a 
card  to  be  the  last  in  the  second  and  third  sets  of  data  cards. 

First  Set 

Card  1:  Column  3  contains  N,  the  number  of  header  cards  that  are 
used  for  comments. 

Next  N  cards:  Comments  in  columns  2  thru  72. 


Second  Set 

Column  1 

Columns  4-6 

Variable 

2 

73 

N7 

2 

85 

NMAX 

2 

96 

NST 

151 

PK(1) 

1 

152 

PK(2) 

Third  Set 

Column  1 

Columns  4-6 

Variable 

6 

PMIN 

9 

UMIN 

10 

EIMIN 

Description 

Number  of  binary  tape  to  be  read 

0.  for  no  tracer  particles  or 
1.  for  tracer  particles 

0.  for  hydrodynamic  computations  or 
1.  for  hydrodynamic  and  strength 
computations 

Problem  number 

Cycle  number  for  restart 

Description 
Minimum  pressure 
Minimum  velocity 

Minimum  specific  internal  energy 
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Column  1 


Columns  4-6 

Variable 

Description 

ii 

RADC 

0.  for  no  radiation  computations  or 
.0137  jerks/cc/kev  for  radiation 
diffusion  computations 

12 

2(12) 

Constant  used  in  flux  limiter  at 
free  surface  in  the  radiation  flow 

16 

THEMIN 

Minimum  temperature  (kev) 

17 

CLITE 

300  cm/shake,  radiation  velocity 

18 

ALPHA 

l#jl  implicit  coefficient  for 
radiation  flow 

21 

DTOPTN 

-1.,  the  code  will  control  the  time 
step  using  the  Courant  and  particle 
velocity  conditions  or 

1.,  the  code  will  use  the  loaded 
constant  time  step 

23 

ESTAB 

Early  time  stability 

24 

EPSLN 

Minimum  flux  (p)  at  free  surface 
when  E  >  EMELT 

25 

FSTAB 

Final  time  stability 

26 

RHOMIN 

Minimum  density  for  cell  to  be 
checked  for  time  change  stability 

27 

BGTHE 

Background  temperature  (kev) 

28 

TCUT 

Use  units  of 

If  0  >  TCUT,  the  stability  (radiation) 
is  checked  if  the  zone  is  either 
cooling  or  heating,  or, 
if  0<TCUT,  the  stability  (radiation) 
is  checked  only  if  the  zone  is  cooling 

37 

S9 

.01,  epsilonics  in  strength 

38 

S10 

Y-component  of  projectile  velocity 
at  t  =  0. 

43 

FEF 

allowable  energy  change,  AE/E,  per 
cycle  for  radiation  flow 

44 

FACT 

Factor  for  increasing  ESTAB  per  cycle 
until  ESTAB  =  FSTAB 
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Column  1 

Columns 

4-6  Variable 

Description 

46 

DMIN 

Problem  will  stop  on  energy  check 
if  | E/E/ cycle |  <  DMIN 

49 

2(49) 

Minimum  DX,  DY,  or  DZ 

Note:  1. 
-1. 

indicates 

indicates 

a  reflective  boundary  condition; 
a  transmittive  boundary  condition. 

2 

59 

LHL 

Left  hydrodynamic  boundary  condition 

2 

60 

LHT 

Top  hydrodynamic  boundary  condition 

2 

61 

LHR 

Right  hydrodynamic  boundary  condition 

2 

62 

LHB 

Bottom  hydrodynamic  boundary  condition 

2 

63 

LRL 

Left  radiation  boundary  condition 

2 

64 

LRT 

Top  radiation  boundary  condition 

2 

65 

LRR 

Right  radiation  boundary  condition 

2 

66 

LRB 

Bottom  radiation  boundary  condition 

2 

67 

NI 

Maximum  number  of  iterations  in  a 
mixed  cell  to  achieve  pressure 
equilibrium 

2 

75 

N9 

0.  for  no  strength  computations  or 

1.  for  strength  computations 

2 

86 

IFRAD 

Maximum  number  of  iterations 
(Newton- Raphson)  for  calculating  the 
temperature  from  the  energy  equation 

2 

96 

NST 

0.  for  no  strength  computations  or 

1.  for  strength  computations 

110 

Z(110) 

1.1*THETA  (background  for  radiation- 
hydrodynamic  computations)  or 

1.1*UMIN  (for  hydrodynamic-strength 
computations) 

2 

116 

LLBB 

Back  hydrodynamic  boundary  conditions 

2 

117 

LHF 

Front  hydrodynamic  boundary  conditions 

2 

118 

LRBK 

Back  radiation  boundary  conditions 
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Column  1 

Columns  4-6 

Variable 

2 

119 

LRF 

2 

123 

11 

2 

124 

12 

2 

125 

J1 

2 

126 

J2 

2 

127 

K1 

2 

128 

K2 

146 

TLP 

147 

DELP 

148 

TDUMP 

149 

DELTD 

150 

TSTOP 

2 

161 

IDEQST(l) 

2 

162 

IDEQST (2) 

2 

163 

IDEQST (3) 

2 

164 

IDEQST (4) 

2 

170 

IDEQST(IO) 

2 

171 

IDKAP(l) 

2 

172 

IDKAP (2) 

2 

173 

IDKAP (3) 

2 

174 

IDKAP (4) 

Description 

Front  radiation  boundary  conditions 

Minimum  active  grid  in  the  X-direction 

Maximum  active  grid  in  the  X-direction 

Minimum  active  grid  in  the  Y-direction 

Maximum  active  grid  in  the  Y-direction 

Minimum  active  grid  in  the  Z-direction 

Maximum  active  grid  in  the  Z-direction 

Time  frequency  for  long  print 

Time  interval  for  long  print 

Time  frequency  for  tape  dump 

Time  interval  for  tape  dump 

Time  for  problem  to  terminate  normally 

Equation  of  state  number  for 
material  1 

Equation  of  state  number  for 
material  2 

Equation  of  state  number  for 
material  3 

Equation  of  state  number 
material  4 

Total  number  of  materials  used  in 
this  problem 

Opacity  routine  for  material  1 

Opacity  routine  for  material  2 

Opacity  routine  for  material  3 

Opacity  routine  for  material  4 
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Column  1 


2 


Columns  4-6  Variable 


Description 


180  IDKAP(IO)  0.  for  no  artificial  viscosity 

computations  or 

1.  for  artificial  viscosity 

computations 

181  OPCON(l)  Linear  viscosity  coefficient 

182  0PC0N(2)  Quadratic  viscosity  coefficient 


183  0PC0N(3)  Maximum  change  in  velocity.  If  the 

change  is  greater  than  0PC0N(3),  the 
artificial  viscosity  coefficient  is 
set  to  0. 


Equation  of  State  Data 


M  =  25 (N  -  1)  where  N  is  the  material  number  (1  thru  4). 

EQSTC(M+1)  thru  EQSTC(M+10)  are  constants  for  Tillotson's  equation 
of  state.  EQSTC (M+ 11)  thru  EQSTC(M+25)  are  constants  that  are  required 
for  any  form  of  an  equation  of  state. 

Column  1  Columns  4-6  Variable  Description 


431 +M 

EQSTC (M+l) 

432+M 

EQSTC(M+2) 

433+M 

EQSTC (M+ 3) 

434+M 

EQSTC (M+4) 

435+M 

EQSTC (M+5) 

436+M 

EQSTC (M+ 6) 

437+M 

EQSTC (M+7) 

438+M 

EQSTC (M+8) 

439+M 

EQSTC (M+9) 

440+M 

EQSTC (M+10) 

441+M 

EQSTC (M+ll) 

442+M 

EQSTC (M+ 12) 

443+M 

EQSTC (M+13) 

P0  (g/cc) 
a 

I0  (jerk/g) 
b 

A  (jerk/cc),  bulk  modulus 
Ig  (jerk/g) 

Ig  (jerk/g) 
a 
6 

B  (jerk/g) 

Not  used 

Not  used 

Not  used 
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Column  1 


Columns  4-6  Variable 


444+M 

445+M 

446*M 

447+M 

448+M 

449+M 

450+M 

451+M 

452+M 

453+M 

454+M 

455+M 


Description 


EQSTC (M+14) 

Not  used 

EQSTC (M+15) 

10"5*P0 

EQSTC (M+16) 

10"6*P0 

EQSTC (M+17) 

IMELT  (jerk/g) 

EQSTC (M+18) 

Yq  (jerk/g) ,  yield  strength  in 

EQSTC (M+19) 

a  in  Y  =  (Y0  +  aP)(l  -  I/IMELT) 

EQSTC (M+20) 

y,  shear  modulus 

EQSTC (M+21) 

p/pQ  failure  criterion 

EQSTC (M+ 22) 

Maximum  pressure  in  Y  function 

EQSTC (M+23) 

Not  used 

EQSTC (M+24) 

Not  used 

EQSTC  Q4+25) 

Specific  heat  (jerk/g/kev) 
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VI.  SUMMARY  AND  CONCLUSIONS 


Several  test  problems  have  been  completed  to  compare  the  numerical 
results  of  an  isothermal  expansion  problem  and  an  impact  problem. 

For  the  impact  problem,  a  plane  strain  (x-y)  impact  (run  with  2DS0IL) 
was  compared  to  a  plane  strain  (x-y)  impact  (run  with  3DS0IL)  with  several 
planes  in  the  z-direction.  These  results  between  two-dimensional  and 
three-dimensional  numerical  techniques  proved  to  be  identical  thereby 
giving  credibility  to  the  2DS0IL  and  3DS0IL  codes. 

A  circular  isothermal  expansion  in  plane  strain  (run  with  2DS0IL) 
was  compared  with  an  equivalent  circular  isothermal  expansion  in  plane 
strain  (run  with  3DS0IL)  with  several  planes  in  the  z-direction.  The 
agreement  was  identical  between  the  two  codes. 

Preliminary  exploration  of  incorporating  an  elastic-plastic  strength 
model  in  the  operator  splitting  mode  indicates  that  it  is  feasible. 

Additional  calculations  will  be  required  to  confirm  the  results. 

The  packaging  concept  (sub-dividing  the  computational  grid  into  N 
problems  and  communicating  the  necessary  information  across  problem 
boundaries)  may  be  necessary  to  enable  one  to  compute  large  three-dimensional 
problems.  Another  possible  solution  would  be  to  pack  several  words  into 
a  single  word.  This  concept  can  double  the  size  of  the  computational  cells 
available  for  any  given  problem. 
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DISTRIBUTION  LIST 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


12  Commander 

Defense  Technical  Info  Center 
ATTN:  DDC-DDA 
Cameron  Station 
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